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ABSTRACT 

We show that energy deposited into an expanding supernova remnant by a highly magnetic (B ~ 
5 x 10 G) neutron star spinning at an initial period of Pj » 2 — 20 ms can substantially brighten 
the light curve. For magnetars with parameters in this range, the rotational energy is released on a 
timescale of days to weeks, which is comparable to the effective diffusion time through the supernova 
remnant. The late time energy injection can then be radiated without suffering overwhelming adiabatic 
expansion losses. The magnetar input also produces a central bubble which sweeps ejecta into an 
internal dense shell, resulting in a prolonged period of nearly constant photospheric velocity in the 
observed spectra. We derive analytic expressions for the light curve rise time and peak luminosity as 
a function of B, Pi and the properties of the supernova ejecta that allow for direct inferences about 
the underlying magnetar in bright supernovae. We perform numerical radiation hydrodynamical 
calculations of a few specific instances and compare the resulting light curves to observed events. 
Magnetar activity is likely to impact more than a few percent of all core collapse supernovae, and 
may naturally explain some of the brightest events ever seen (e.g., SN 2005ap and SN 2008es) at 
L > 10 44 ergs s~L 



Subject headings: radiative transfer - stars: neutron - supernovae: general 
(SN 2005ap, SN 2008es, SN 2007bi) 



supernovae: individual 



1. INTRODUCTION 

Studies of soft gamma-ray repeaters and anoma- 
lous X-ray puls ars reveal that ~ 10% of newly born 
neutron stars (jKouveliotou et all [1998) have dipole 
magnetic fields as high as B ~ 10 14 — 10 15 G 
for more than 1000 ye ars after their birth (see 
iWoods fc Thompson! r2006). These "magnetars" rotate 
at periods of P = 5 - 12 s at an age of 1000 - 10,000 
years. Such highly magn etized neutron stars (NSs) 
were theoretically predicted ([Duncan fe Thomp son 1992; 
iThompson fe Duncan! Il993h , and most of their activity 
(both sporadic and persistent) must be powered by the 
decay of these large magnetic fields. 

What remains unknown is just how highly magnetized 
and rapidl y rotating these magnetars may be at "birth" 
Many (see iBodenheimer fe Ostrik"er] [l974: Wheeler et alj 
l2000t IThompson et al.ll2004h have investigated the pos- 
sible impact on the central engine when the magnetar is 
so rapidly rotating (1-3 ms) and magnetized that its sub- 
sequent spin-down can power the explosion. Cases this 
extre me may also be sources for ultra-high energy cosmic 
rays (jArons 2003T) or deposit enough energy in the col- 
lapsing stejla£_ejiyek2p_e^^ 

terior (jUzdenskv fc Ma cFadvcn 2007]; IBucciantini et all 
2009) for the production of a collimated relativistic flow 
needed for gamma-ray bursts. Such events depend on 
the combination of rapid rotation and high B to achieve 
a measurable effect during the few seconds critical to the 
core collapse mechanism. 
Building on the work of IGaffetl CEIZZlIB), 
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found that weaker magnetic fields and less extreme spins 
(that do not alter the explosion mechanism) can dramati- 
cally impact supernovae light curves, competing with the 
decay of radioacti ve 56 Ni and thermal energy in the ex- 
panding envelope. Macda ct al. (20071) previou sly raised 
this p ossibility for the Type lb SN 2005bf, and lWooslevI 
(2009) has independently shown their relevance as well. 

We show in fJ21 that when the timescale of the mag- 
netar spindown, t v , is comparable to the effective ra- 
diative diffusion time, td, the resulting peak luminos- 
ity is Lp 0a k ~ Eptp/t^, where E p is the magnetar rota- 
tional energy. Magnetars with 10 13 G < B < 10 16 G 
and Pi = 1 — 30 ms can produce L pca k > 10 42 erg s~L 
We discuss the dynamics of the energy injection in fJ3] 
and show that the magnetar blows a central bubble in 
the SN ejecta, forming a dense inner shell of swept- 
up material which affects the spectroscopic evolution. 
In 21 we derive analytic expressions for the luminos- 
ity, Xpeak, and duration, i pC ak, of magnetar powered su- 
pernovae. We confirm these formulae with numerical 
radiation-hydrodynamical calculations, and show how 
they can be inverted to infer B and Pi from a given light 
curve. We close in *j5]by discussing observed core-collapse 
SNe that may be p owered this way, esp ecially the ultra- 
bright SN 2005ap (fQiumbv et al.N2007l) and SN 2008es 
(|Gezari et al.ll2009t iMTller et al.ll2009D . 

2. MAGNETAR HEATING: SIMPLE ESTIMATES 

In the simplest model, the core collapse mechanism 
has ejected an envelope of mass M c j at a velocity i>t 
from a star of initial radius Rq. Within a few ex- 
pansion times, t G ~ R /v t , this ejecta will be under- 
going self-similar adiabatic expansion, with an internal 
energy E int ~ E sn (R a /R), where E sn ~ M cj v%/2 and 
R ~ v t t is the remnant size. In the absence of mag- 
netar (or 56 Ni) heating, adiabatic expansion continues 



2 



Kasen & Bildsten 



until the remnant is as old as the effective diffusion time 
td ~ (nMcj/vtc) 1 / 2 , where k is the opacity, after which 
the entropy is lost. Such thermally powered light-curves 
(e.g. Type Hp's) have a luminosity L t h ~ E sn t c /t\. The 
large amount of adiabatic expansion that has occurred 
by the time t ~ td leads to low luminosities. 

Now consider the impact of late time (t 3> t c ) en- 
ergy injection from a young magnetar with radius R ns = 
10 km and initial spin fij = 27r /Pi. The magnetar rota- 
tional energy is 



E„ 



2 x 10 50 P^ 2 ergs, 



(1) 



where P 10 = Pi/ 10 ms and we set the NS moment of in- 
ertia to be I ns = 10 45 g cm 2 . This magnetar loses rota- 
tional energy at the rate set by magnetic dipole radiation 
(with the angle, a, between rotation and magnetic dipole 
fixed at sin 2 a = 1/2), injecting most of the energy into 
the expanding remnant on the spin-down timescale 



tr 



6/nsC 3 



= 1.3B^ 2 P? yr, 



(2) 



where P14 = P/10 14 G. To input this energy at a time 
^5 requires a minimum B field of 

B > 1.8 x 10 14 P10 k" 1/4 Af 6 ~ 3/8 ^( 8 G, (3) 

where k cs = k/0.2 cm 2 g _1 , M 5 = Af cj /5 M Q and E 51 = 
P sn /10 51 ergs -1 . The required fields are in the magnetar 
range. This late time entropy injection resets the interior 
energy scale to Ei nt ~ E p and overwhelms the initial 
thermal energy when E p > E sa (t e /t p ). Thus even low 
magnetar energies E p < E sn play an important role. The 
resulting peak luminosity is 



^pcak 



Ept p 
t 2 



5 x 10 43 B^ 2 K - 1 M 5 3/2 El( 2 erg s 



( 4 ) 

which is primarily a function of the magnetic field. This 
shows that L pcak - 10 43 - 10 45 erg s -1 SNe can be 
achieved from magnetars with B^ = 1 — 10 and initial 
spins in the P, = 2 — 20 ms range. 

3. HYDRO DYNAMICAL IMPACT 

Our simple estimate ignores the details of how the de- 
posited energy is distributed throughout the interior of 
the expanding SNe remnant. Since the dissipation mech- 
anism for the pulsar wind in this medium is poorly under- 
stood, we assume the injected magnetar energy is ther- 
malized spherically at the base of the supernova ejecta. 
The remnant is assumed to be in homologous expansion 
with a shallow power law density structure in the interior 
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where vt — (2E sn /M e j) ' is the characteristic ejecta ve- 
locity, and the density falls off sharply above vt ■ 

The central overpressure caused by the energy deposi- 
tion blows a bubble in the SN remnant, similar to the 
dynamics studied in the context of pulsar wind neb - 
ulae (e.g.. lChevaiienll977t [Chevalier fc Franssonlll992fl . 
As this bubble expands, it sweeps up ejecta into a thin 




0.4 0.6 0.8 1.0 
velocity coordinate (10 9 cm/s) 

Fig. 1. — Radiation-hydrodynamical calculations of the density 
(top) and temperature (bottom) of magnetar energized supernovae, 
one month after the explosion. The supernova had M e j = 5 Mq 
and E Bn = 10 51 ergs. The magnetar had t p = 10 sec and various 
values of E p , labeled in units of 10 ergs. The dashed line in 
the top panel shows the unperturbed density structure, taken from 
Equation JSJ. 

shell near the leading shock, leaving the hot, low den- 
sity interior evident in the TD radiation hydrodynam- 
ical calculations of Figure [TJ In multi-dimensional cal- 
culations of pulsar wind nebulae, Rayleigh- Taylor insta- 
bilities broaden the shell and m ix the swept-up material 
(|Junlll998l : iBlondin et aT1l2001h . 

The bubble expansion will freeze out in Lagrangian 
coordinates when the leading shock velocity becomes 
comparable to the local velocity of the expanding SN 
ejecta. The postshock pressure is P = 2jp v 2 /(l +7) = 
(8/7)po^ 2 for a strong shock, and the pressure of the 
energized cavity is P rs E p /3V, where V is the volume, 
implying a shock velocity v 2 = 7E p /32irR 3 po. The shock 
becomes weak when v s w R/t, which determines the final 
velocity coordinate of the dense shell 
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16(3 - 5) E s 



1/(5-5) 



for E p < E s -_ 



The weak dependence on E p , v sh oc Pp /4 , for 6=1, 
places v s h near i> t . The total mass swept up in the shell 
isM sh =M cj (« t /D sh ) 3 ^. 

The magnetar does not affect the dynamics of the outer 
layers of the SN ejecta unless E p > P sn , in which case the 
bubble expands beyond v t and accelerates more rapidly 
down the steep outer density gradient. Essentially all of 
the ejecta is then swept up into the shell at a final shell 
velocity 



(6) 



Vsh 



v t [l + E p /E sn ] 



1/2 



for E p > E s 



(7) 



Both estimates for w s h assume no radiative losses. 

The presence of a dense shell has consequences for the 
supernova spectra. Initially the photospheric velocity, 
w p h, as measured from the Doppler shift of absorption 
line minima, decreases with time as the outer layers of 
ejecta become transparent. Once v p h has receded to the 
shell velocity, however, it will remain constant (Figure [2l 
bottom panel). The spectra will then be characterized 
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Fig. 2. — Radiation hydrodynamical calculations of magnetar en- 
ergized supernovae with M c j = 5 Mq, E sn = 10 51 ergs, and a den- 
sity structure given by equation J5]l. The magnetar had Pi = 5 ms, 
and various magnetic field strengths. Top panel: Bolometric light 
curves. The dashed line shows, for comparison, the energy depo- 
sition from 1 Mq of 56 Ni. Middle panel: Effective temperature. 
Bottom panel: Velocity of the e~ scattering photosphere at r = 1. 



by relatively narrow but blueshifted absorption features, 
and the spectral evolution will be notably slow. The shell 
becomes optically thin to electron scattering at a time 



t T=1 = 326 M t 



£T :/ V/ 2 



days. 



(8) 



Recombination may hasten this transition. The photo- 
spheric velocity drops suddenly to zero after t T —i, how- 
ever because the line opacity in the shell remains opti- 
cally thick longer, the spectral features may not display 
any dramatic change for some time after. 

4. LIGHT CURVES 

We now derive analytic expressions for the peak lu- 
minosity of a magnetar powered SNe using a one-zone 
model for the whole remnant. The internal energy, £j n t, 



is governed by the first law of thermodynamics 



dt 



-P 



dV 
~dt 



+ Lp(i) - L e (t), 



(9) 



where L p is the magnetar luminosity and L e the radi- 
ated luminosity. We assume that the magnetar energy is 
thermalized throughout the remnant, and that radiation 
pressure dominates, P = Ei nt /3V. When the volume 
increases as V oc i 3 , equation (|9|) becomes 



1 d 



[E int t] = L p (t) - L e {t). 



(10) 



t dV 

The radiated luminosity, L c , is approximated from the 
diffusion equation 

L c _ c 8E int /V c E int /V 



AttR 2 Sup dr 



3np R 



(11) 



and rewritten using R = Vft, defining the effective diffu- 
sion time, id 



L c = 



E\„tt 



where — 



3 M ci K 

47T VfC 



1/2 



(12) 



where we take Vf = [(E p + E sn )/2M e j] 1 / 2 as the final 
characteristic ejecta velocity. For the simple case where 
the magnetar injects a constant luminosity L p = E p /t p 
over a time i p , and then shuts off, we find 



L e (t) 



L c (t) = 

t p 



-* 2 /2*dl 



* 2 /2*d [ e *p/2*d 



t < u 



t > u 



(13) 



This light curve peaks at a time t pi then declines on 
the characteristic time scale id- For t p <g; id, -^pcak = 
-E p ip/2ij, similar to the estimate in §2. When i p 3> id, 
we find L pcak = E p /t p . 

More generally, the energy input from the magnetar 
persists for i > i p , and is given by the spin-down formula 



L P (t) = 



EL 



I - 1 



i p (1 + i/ip)' 



(14) 



where / = 2 for magnetic dipole spin down. The energy 
input at late times may not be dynamically important, 
but it enhances the luminosity by continually heating 
the ejecta in a manner similar to the decay of 56 Ni. No 
simple analytic solution for the light curve exists for the 
general form of L p (t), but since radiative losses are mini- 
mal for times i < id we can derive approximate relations 
by solving equation (fTOf for the case L e = 0. The re- 
sulting internal energy can be evaluated at time id in 
equation fj 1 21) to estimate the peak luminosity 



L 



peak 



L 



peak 



r E ptp 
t 2 

j.E p t p 



111 



1 + r 



id + ip 



i 2 

'd 



1 - 



id/t P (Z-l) + l 



(1 + td/tp) 1 - 1 



1 = 2, 

I > 2, 
(15) 

where the correction factor / will be calibrated by com- 
parison to numerical simulation. In general, L pea k de- 
creases as I increases, as more of the energy is deposited 
at earlier times and suffers greater adiabatic losses. 
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Fig. 3. — Required B and Pi needed to achieve a given £ pea k- The 
lines are contours of constant -t/ pea k assuming E an = 10 51 ergs and 
M e j = 5 Mq (solid) or M ej = 20 M Q (dashed) from equation lT5t . 
Regions to the right of the knee have t p < t^, whereas regions to 
the left of the knee have t p > tj. The horizontal dotted line shows 



where E v 



10° 



At the peak of the light curve, the radiated lumi- 
nosity equals the instantaneous magnetar luminosity, 
-kpeak = Lp(t pea k), the general expression of "Arnett's 
law" (Arnett 1979). This follows from equation (fT0|) . 
since equation (|12[) implies that the maximum of L e oc- 
curs when d{E ln tt) / dt = 0, yielding the time of maxi- 
mum in the light curve 



t 



peak 



= U 



(l - 1)E P 



- 1 



(16) 



^pcak^p 

For tp <C td the light curve peaks at i pea k « 
td/ _1 ^ 2 [ln(td/i P ) — l] -1 / 2 (assuming / = 2), whereas for 

t p 3> id the peak occurs at i pca k ~ tp(y/Wf ~ -0- 

Figure [2] shows 1-D radiation hydrodynamical calcu- 
lations for M e j = 5 Mq,E si1 — 10 51 erg, and central 
magnetars (/ = 2) with Pj = 5 ms. A grey opacity 
k = 0.2 cm 2 g _1 was assumed. The simple one zone 
model works remarkably well at predicting L poak and 
ipoak and comparison with the numerical models fixes 
the value of / = (/+ l)/2. At late times (t > t T= i) when 
the SN becomes optically thin, the light curve tracks the 
magnetar luminosity, L ~ i~ 2 , which is similar to the 
curve of 56 Co decay. Late time measurements of the 
bolometric light curve could discriminate the two energy 
sources, though it is not clear that the assumptions of 
complete thermalization and constant I = 2 spindown 
will hold at these late times. 

In Figure [3l we use equation (| 1 5() to find the locus in 
the Pi — B space (assuming I — 2) needed to reach a 
certain -L pca k in a supernova with E sn = 10 51 ergs and 
M e j = 5 M Q or M ej = 20 M . A larger M oj increases 
id, which reduces L pea k for a given set of magnetar pa- 



FlG. 4. — The dependence of L pca k an d ipeak on the initial 
magnetar spin and B field. The solid lines are for fixed B14 = 
100,30,10,3,1,0.3 and 0.1 and varying spin period, whereas the 
dashed lines are for a fixed Pi = 1, 3, 10 and 30 ms and varying B. 
This calculation assumed E Bn = 10 51 ergs and M e j = 5 Mq. 

rameters. Magnetars with Pi < 5 ms (below the dotted 
line) dump enough energy to increase the ejecta veloc- 
ity, shortening t^- The lines merge for low B as they 
asymptote to L pca k — > L p . 

We can also invert the problem and use the measured 
values of L pea ^ and £ pea k for an individual supernova to 
infer B and Pi. Figure 2] uses equations (TT5)) and (|16[) to 
illustrate how L pe ak and t pea k vary with B and Pj. This 
"mapping" allows for an assessment to be made of the 
magnetar 's properties and illuminates which numerical 
calculations should be done. We placed the observed 
values for 2008es on this plot, motivating the numerical 
results we show in the following section. A different plot 
would need to be made for different M e j and E sn . 

5. DISCUSSION AND CONCLUSION 

We have shown that rotational energy deposition from 
magnetar spin-down with initial spin periods < 30 ms 
can substantially modify the thermal evolution of an ex- 
panding SNe remnant. For magnetars in this range, the 
peak luminosity reaches 10 42 — 10 45 erg s _1 (A/boI = 
— 16.3 to —23.8), substantially impacting the typical 
core-collapse SNe lightcurve, whether it is a Type II or a 
Ibc event. The highest luminosities occur when t p ~ t<i, 
in which case the total energy radiated in the light curve 
is i? ra d ~ L pca k^d ~ Ep/3. The maximal spin of a NS is 
around 1 ms, so E Ta d cannot exceed ~ 10 52 ergs; super- 
nova radiating larger energies can not be explained by 
this mechanism. Though we know that ~ 10% of core 
collapse events make magnetars, we do not know the dis- 
tribution of initial spin periods, so the prevalence of light 
curve dominance is difficult to predict. 

For stars with remaining hydrogen, magnetar injec- 
tion may explain the brighter (Mb ~ —19) subclass 
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Fig. 5. — Bolometric light curve calculations of magnetar ener- 
gized supernovae compared to observed events. A constant opacity 
k = 0.2 g cm -2 is assumed. Black circle s show V-band ob serva- 
tions of the luminous Type IIL SN2008es HGezari et aLll2009l> with 
an assumed rise time of 25 days. Red squares show R-ba nd obser- 
vations of the Type Ic SN 2007bi UGal-Yam et al.ll2009h with an 
assumed rise time of 50 days. 



of Type II-L SNe noted by iRichardson eTail (|200l . 
i.e. 1961F, 1979C, 1980K, 1985L. The light curves of 
these events are difficult to explain in standard explosion 
models unless extreme progenitor radii (R > 2000 Rq) 
are assumed (Blinniko v k. Bartunovl 119931 ). Figure \5\ 
shows that a magnetar with relatively modest rotation, 
Pi = 10 ms, in a M C j = 5 Mq supernova can reach 
similar luminosities. Events brighter than Mb i = —21 
(L > 8 x 10 43 er gs" 1 ), such as the ultrabright Type 
II-L SN 2005ap (iQuimbv et al l 12001 and SN 2008es 
(jGezari et al.ll2009t IMiller et al.ll20'09f ) require an initial 
magnetar spin of < 5 ms. Motivated by Figure [U we 
found an excellent fit to the SN 2008es light curve with 
B14 = 2, Pi = 2 ms and M e - 3 = 5Mq. Such rapidly ro - 
tating magnetars must be rare, as lVink fc Kuiperi (j2006l ) 
found that the galactic supernovae remnants of known 
magnetars were explained with typical explosion energies 
of 10 51 ergs. This rarity is consistent with the specific 
volume rate of these ev ents; current estimates put them 
at no more than ~ 1% (|Miller et al.ll2009l : IQuimbv et all 
2009) of the local core collapse rate. 



Debate remains ( seelKlose et al.l 12004 IGaensler et al.l 
120(51 iDavles et all 12001 " as to whether magnetars are 
preferentially formed from the most massive stars that 
collapse to NSs. If so, then we might see a prevalence 
of magnetar dominated light curves amongst the Ib/c 
SNe, which may partially explain the wide light curve 
diversity note d in this class. Some e xtreme SN Ic, such 
as SN2007bi (IGal-Yam et all I2TM lYoung et all [200§) 
and SN 1999as (jKnop et al.l ll999 s |. which remained very 
bright for a long time, have been claimed to be the 
pair instability explosion of a ~ 100 Mq star produc- 
ing nearly 5 Mq of 56 Ni. Figure [5] shows that the 
light curve could alternatively be explained for a su- 
pernova with M e j = 20 Mq forming a magnetar with 
£14 = 2, Pi = 2.5 ms. The spectra of SN 1999as 
also revealed a slowly evolving photospheric velocity and 
narrow, blueshifted absorption fea tures, sugges tive of a 
dense shell like that predicted here (|Kasenll2004f ). On the 
other hand, the magnetar model may have trouble repro- 
ducing the strong iron emission lines seen in the nebular 
phase spectrum of SN 2007bi. 

Our initial investigations have revealed that if an 
appreciable fraction of highly magnetic NSs are born 
rapidly rotating, then we should find evidence for them 
in the plethora of super novae surveys, su ch as the Palo- 
mar Transient Factory (|Law et~al][2009l ). Many open 
questions remain on the theoretical side, especially how 
the outgoing pulsar wind thermalizes in the remnant, 
whether there are substantial Raylcigh- Taylor instabil- 
ities, and how these could manifest themselves in the 
observed spectra both at late times and during the pho- 
tospheric phase. Our work has outlined the regimes of 
relevance, and will guide future large scale computations 
through parameter space in an informed manner. 
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